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Summary A hybrid discrete-continuous model of layered superconduc- 
tors with interlayer Josephson couplings of arbitrary range is constructed. 
The conditions required by gauge invariance and thermodynamical stability 
of the model are determined. Some important special cases, in particular 
transition to the classic Lawrence- Doniach model, are discussed. The con- 
ditions for presence of alternating solutions and the posibilities to describe 
such states of a superconductor by the continuum or bi-continuum models 
are examined. The enhancement of superconductivity caused by the presence 
of higher order Josephson couplings is shown. 

1 Introduction 

Most of the high-temperature superconductors like e.g. YBCO or BSCCO 
have a layered structure. Such a strongly anisotropic situation results in 
the fact that the material properties and behaviour of the fields in direction 
(say z axis) orthogonal to the layers is totally different from the behaviour 
in directions paralel to layers. [4, 1, 8]. The idea originally proposed by 
Lawrence and Doniach (LD) [7] is to consider the Ginzburg-Landau order 
parameter ip as a function of two continuous variables (say x and y) and 
one discrete variable n - the index of the layer. The form of the free-energy 
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functional proposed in [7] and, after a modification, presented e.g. in [5] takes 
into account the Josephson coupling between the nearest neighbour layers. 
The exact solutions for this case are given in [6]. The higher grade hybrid 
model, proposed in this paper, admits also couplings between more distant 
neighbours - up to a (given, but arbirary) range K. 

Let us look at the layered superconductor as a one-dimensional chain of 
atomic planes with Josephson's bonds between them. Such bonds will be 
called J-links. The interplanar distance will be denoted by s. We assume the 
following convention for indexing the planes and links. If we locate the point 
z = at an atomic plane, then the z-coordinate of any plane, equal ns, may 
be represented by the integer n, while the z-coordinate of the center of any 
interplanar gap, equal Is, by the half-integer /. Choosing the point z = at 
the center of an interplanar gap - we index the planes by half-integers and 
the gaps by integers. 

2 The hybrid model of grade K 

Let us consider the free-energy functional T for a layered superconductor. 
We shall denote by ip n the order parameter associated to the layer indexed 
by the number n. Its complex conjugate (c.c.) will be denoted by ip n . The 
symbol m ab and m c will denote the in-plane and tunneling effective mass of 
superconducting current carriers, respectively. We start from the free energy 
functional of the following form 

F = F + F s + ±-j B 2 d 3 x. (1) 

The term JF describes the normal state, while T s the superconducting one. 
The supercoducting term is composed of two parts: 

T s = F v + Tj, (2) 

where the part 

T p = Y,Fn (3) 

n 

describes the contribution of atomic planes, while Tj corresponds to inter- 
planar Josephson's bonds. For any plane indexed by n the free energy F n 
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has the 2D Ginzburg-Landau form (in general, the parameters can depend 
on n) 

F n = sJ dxdy{^-\{^) n \ 2 + aolVnl 2 + ^lV>n| 4 }, (4) 

where we have introduced the 2-dimensional continuous operator D (covari- 
ant derivative) 

ie* 

D p = d p - foAp, p = x,y. (5) 

The form (4) of the functional T v already ensures its invariance with 
respect to the gauge transformation 

A ^A' = A + VA, ^ n ^ = ^e^ A , ^ n - < = ^e"^ A (6) 

(2-dimensional A and V for this case). The standard variational treatment 
of the functional T s with respect to A p , p — x,y, gives the standard in-plane 
components of the supercoducting current 

ie*k - - e* 2 

jp = ~7, {^ndplpn - Ipndp^n) A>|^n| , p = X,y. (7) 

2m ab p m ab c 

Let us now construct the term Tj in the free-energy. In general it is a 
functional which can depend on all V'ruV'ra an d on the vector potential A. 
Consider first the global gauge transformation, i.e. A = const in (6). The 
invariance condition for Tj implies that its density Fj fulfills the relation 

dFj i dFj i n <*\ 

77— - 77— = 8 

Olpn dl{j n 

for each n separately. That means that the functional Tj depends on the 
fields ip only through the combinations ip n ipk- 

Now we shall consider the local gauge transformation with A depending 
only on the variable z, but first let us introduce the following new variables 

i>n = ^ n e ipin (9) 
(and the appropriate complex conjugate), with 

g* r(n+l)s 

p ln = — A z dz. (10) 

TlC Jns 



The vanishing of the variation of Tj with respect to A gives the condition 

which, together with the condition of global gauge invariance, implies explicit 
independence of Tj of A z . 

The simplest gauge invariant expression for the energy of J-link between 
n-th and (n + g)-th planes will have the following form 

€ qn = ^{Cqni'n^n + Vqn4>n+q^n+q - {iqn^n^n+q^^ + C.C.)}, (12) 

where the exponent p qn is defined by the formula 

g* r(n+q)s 

Pqn = T~~ / A g dz. (13) 

tlC J ns 

Let us note that for the model invariant with respect to the time reversal we 
have 7 gn = In general, the coupling parameters (, i], 7 as well as in-plane 
parameters a and (3 can be different for different planes (which is the case 
for superconductors composed of various atomic planes). However, in this 
paper we shall consider the array of identical planes, hence the parameters 
will not depend on index n. That implies that i] q = ( q . Hence instead of (12) 
we shall use 

e 9 n = \{( q (\i>n\ 2 + \^n+q\ 2 ) ~ Iq^n+qe^ + C.C.)}. (14) 

Let P denote the (finite or infinite, but ordered) set of all indices of planes, 
and let Q = {1,2,..., K}. The planes connected with the n-plane by Joseph- 
son coupling select the following subset of Q 

P n = {qeQ: (n + q)eP}. (15) 

Thus, the gauge invariant functional Tj for the hybrid model of grade K has 
the density of the following form 

Fj = E E (16) 

neP qeP n 

with e qn given by (14). The coupling parameters ( q and 7 9 vanish for q > K. 
Every J-link is represented in (16) by exactly one term. 
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3 Comparison with anisotropic GL model 



To compare our hybrid model (HM) with the continuum GL model, we shall 
consider the infinite medium. In that case the summation index may be 
shifted by the integer q. Moreover P n = Q. This implies that 

f j = \ E E{2(c - i q )m 2 + 7 ff i - ^i 2 } ( 17 ) 



n q 



For very weak field A z and very small dependence of ip n on n we have the cor- 
respondence rules which allows us to pass from the hybride to the continuum 
models. 



qs 



(18) 



l-ifcqsA z . 

Applying the rules to the functional (2) with (4) and (17) one obtains 



J d 3 x{ 



ft 



2m, 



ab 



|DVf + c*o|V| 2 + o/W + 



1 



+ E[(C,-7,)l^| 2 + ^V 7g | J D 3 ^| 2 ]}. 



(19) 



Hence, we have the following relation between m c - the effective mass in 
z— direction (in anisotropic GL model) and the coupling parameters 7 g : 

k = wV<- (20) 

The presence of J-links modifies also the parameter a to the form: 

a = a + J2((q-%)- (21) 



4 Field equations 

By computing the variation of the functional T with respect to A z one obtains 
the Maxwell equation for the z— components of current density and curl if 

-J( z ) = ^(cm\H) z , (22) 
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where 

se* - 

J ( Z ) =^EE {l q ll>ni>n +q e- Wqn - C.C.} Xg n(z), (23) 
mC neP qeP n 

P n is given by (25), and the quantity p qn by (13). The symbol Xqn( z ) denotes 
the characteristic function of the interval [ns, (n + q)s]. Let us note that 
for any layer ns < z < (n + l)s the exppression J(z) does not depend on 
the value z; only the ends of the interval are important. To better see the 
structure, let us first extend the set Q on the negative values 

Q = {-#,.., -2, -1,1, 2, ..,#}, (24) 
and introduce the symbols 

„ ._ / 1 ^ (n + q) eP , . , 

qn \ otherwise . { > 

Then the expression for Josephson current J; describing tunneling across the 
interplanar gap indexed by I (half integer if P contains integers) will have 
the form 

se* - 

J l = ^ E E {lq^qn^n+ q e- Wqn ~ C.C.}, (26) 
ltlC qeQ neP lq 

where 

P iq = {ne P : n < I < n + q}. (27) 

By computing the variation of the functional T s with respect to ip n , one 
obtains the equations 



-^D^n + + P\lpn\ 2 1pn+ 

(28) 

-\E q eQlq(^qn^n+ q e iPqn + O -qu^u-q^^*) = 0, 

where instead of ceo we have introduced 

a = a + E <VC« ( 29 ) 

depending of n, for finite P. 
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5 The ground states 



Let us now consider the plane-uniform states of HM in the absence of mag- 
netic field. The order parameter is then independent of the in-plane variables 
and the net supercurrents vanish. In detailed calculations we shall focus our 
attention on the grade K = 2, which seems to be sufficiently illustrative to 
grasp the idea on what is going on. The generalization to grades K > 2, 
although more complicated algebraically, is straightforward. In the region 
far from the boundary all the coefficients a qn = 1. For K = 2 the condition 
of vanishing Josephson current is equivalent to 

7i(^„^„+i - c.c.) + 7 2 (^„^„+2 + ^n-iVWi - c.c.) = 0. (30) 
and the equations (28) take the form 

aipn + P\iJ n \ 2 ^n - ^[7l(VWl + VV-l) + 72(VW2 + i>n-2)] = 0. (31) 

We shall look for solutions with constant amplitude and diference of phase 
between nieghbouring atomic planes, so we use the ansatz 

= Ce me . (32) 

The result is the equation 

a + - 71 cos -72 cos 20 = 0, (33) 

with the condition 

7! sin + 27 2 sin 20 = 0. (34) 

Solving (34) with respect to 0, we obtain 3 variants: (al) = 0, (a2) = it, 
and (a3) cos0 = — (71//472). The solution C to (33) has the form C 2 = 
— a*// (3, where a* depends on the variant. The variant (al) implies the 
uniform solution to (31): 

= C, (35) 

with 

a* = a + Ci + C2 - 7i - 72, (36) 

what is the case of equation (21) for grade K=2. In the variant (a2) the 
solution to (31) has the alternating form 

= ±C, (37) 
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with 

a* = a + Ci + C2 + 7i - 72- (38) 

Finally, in the variant (a3), the solution exists if the parameters 71 and 72 
fulfill the relation 

| 7l |<4| 72 |. (39) 

Then the parameter a* is connected with the coupling constants by the 
formula: 

a* = a + Ci + (2 + 72(1 + ^2). (40) 

872 

There are two independent solutions 

^j n = Ce ±d , 9 = arccos(-^). (41) 

72 

They will be referred to as the phase modulated states. The solutions de- 
generate at the extremities |7i| = 4| 72 1. 

So far we confined our discussion to the existence of solutions which could 
serve as candidates for the ground state. The question of stability of the 
solutions will be addressed in the next section. 

Let us note that the condition sin# = admits the solutions (35) and 
(37) for hybrid model of any grade K. Contrary to that, the analogues of 
the conditions (a3) and (39) can deliver, depending on the grade and the 
coupling parameters, any number from to 2{K — 1) modulated solutions. 



6 Stability 

To examine the stability of the solutions presented in the previous section, 
we shall analyze the Hessian matrix of the free energy T s describing small 
deviations from the ground state. The problem reduces to examining the 
function 

£(C,0) =aC 2 + -[3C 4 -C 2 ( 7l cos 9 + 72 cos 29). (42) 
2 

The solutions found in the previous section are stationary points of £. If 
72 = 0, then the stability of the solutions depends on the sign of 71. If 
72 7^ then the sign of the respective eigenvalue depends on the values of the 
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parameters 71 and 72. For 6 different from and tt the dependence has form 
according to the function 

/(7i,72) = 72(71 + 472) (71 - 472). (43) 

The straight lines 71 + 472 = and 71 — 472 = divide the plane 71, 72 into 
4 regions (see Figure 1). 

As explained above, both the uniform and the alternating solutions always 
exist. However, in the region 

(A) : 71 > 0, -71 < 4 72 < 71, (44) 

only the uniform solution (35) is stable. 




Fig. 1. The regions of stability: uniform (=) , 
alternate ± , and phase modulated (~) , solutions 



On the other hand, in the region 

(B) : 72 > 0, -4 72 < 71 < 472, (45) 

both the solutions (35) and (37) are stable (the uniform one in the right and 
the alternating one in the left part of the region), while in the region 

(C) : 71 < 0, -71 < 4 72 < 71, (46) 
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only the alternating solution (37) is stable. The region 



(D) : 72 < 0, 4 72 < 7i < -472, 



(47) 



excludes the stability of both the uniform and the alternating solutions but, 
in contrast to that, ensures the existence and stability of the modulated 
solutions (41). In the region (45) the modulated solutions exist but are 
unstable. 

Let us note that in the regions (45) and (46) of stability of the alternat- 
ing solution (37) one can apply the construction of bi-continuum solution 
presented in [9]. 

7 Enhancement of the superconductivity 

The asociation of the formulae (36), (38) and (40) with the regions of stability 
of the ground state shows that, for suitable relations between the coupling 
constants 7 i and 72, one can make the parameter a* more negative than 
a . In consequence, the 3D superconductivity can turn out enhanced with 
respect to the 2D superconductivity in the atomic planes. Such a possibility 
has been indicated in literature [2]. In further discusion we shall count the 
enhancement with respect to the origin located at a given by (29). It is 
convenient to introduce the polar coordinates 7 and ip in the plane 71, 72: 



71 = 7 cos (p, 72 = 7 sin (p, 



(48) 



and the notation 



ipo = arctan(-). 



(49) 
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Fig. 2. Enhancement of superconductivity 
by Josephson currents: Aoc vs. cp 

The quantity Aa = «*-«asa function of the coupling angle p is plotted 
in Figure 2 (the numerical values are computed for 7 = 1). 
In the uniform state (36) we have 

3 

Aa = ^7 sin((/9 — -7r), — (p <<p <n — ip Q . (50) 

The minimum Aa (hence the maximum enhancement) is reached at p = 7r//4 
and equals —y/^j. 

In the alternating state (38), in turn, one obtains 

3 

Aa = V2^y sin(ip + -tt) , +ip <ip < tt + ip (51) 

with the minimum value — v / 2 _ 7 reached at ip — jir. 

The enhancement for the phase modulated state (40) is represented by 

Aa = 7siny?(l + -ctgV), -n + V?o < V < - ¥o- (52) 
8 

In this case the minimum equals —7 and is reached at ip — —it/ /2. Hence, 
the maximum enhancement in this case is smaller than in the uniform and 
alternating states. 

8 Special cases 

The hybrid model of grade (HM with K=0) reduces formally to uncoupled 
(2 dimensional) Ginzburg-Landau model for each atomic plane. The HM 
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with K=l has in general two coupling parameters ( and 7. If they are equal 
to one another, one obtains the Lawrence-Doniach model with parameter 
a = a , and the effective mass in z-direction given by the formula 

1 s 2 

— = |*7- (53) 

Due to 7 > the LD model has only one stable ground state solution, namely 
the uniform one. The interplanar coupling gives neither enhancement nor 
suppression of the critical temperature. Although the Josephson current 
coupling places the model on the enhancement side, the efect is precisely 
annuled by ( = 7. 

Let us note that, in the enhancement mechanism discussed above, the 
2D superconductivity of the planes is not a prerequsite for the 3D supercon- 
ductivity of the array of the planes. In fact, one can obtain the negative 
a* starting from positive a - This is in concordance with ideas expressed in 
Anderson's discussion of his Dogma V in [3]. 

9 Acknowledgements 

This work was supported by the Science Research Committee (Poland) under 
grant No. 5 T07A 040 22. 



References 

[1] A. A. Abrikosov, Fundamentals of the Theory of Metals, North-Holland 
1988 

[2] P. W. Anderson, Intralayer tunneling mechanism for high-T c super- 
conductivity: comparison with c-axis infrared experiments included in: 
P.W.Anderson, The theory of superconductivity in the high-T c cuprates, 
Princeton University Press 1997 

[3] P. W. Anderson, The theory of superconductivity in the high-T c cuprates, 
Princeton University Press 1997 



12 



[4] G. Blatter, M. V. Feigel'man, V. B. Geshkenbein and A. I. Larkin and 
V.M. Vinokur, Vortices in high-temperature superconductors, Reviews 
of Modern Physics 66, 1125(1994) 

[5] J. B. Ketterson and S.N. Song, Superconductivity, Cambridge Univ. 
Press 1999 

[6] S.V. Kuplevakhsky, "Vortices and the Meissner effect in stacked junc- 
tions and layered superconductors: Exact analytical result, Phys. Rev. B 
63, 054508 (2001) 

[7] W. E. Lawrence and S. Doniach, Theory of layer structure superconduc- 
tors, in Proceedings of the Twelfth International Conference on Low 
Temperature Physics, Kyoto, ed. E.Kanda, Academic Press of Japan 
1971, p.361 

[8] D. Rogula, Dynamics of magnetic flux lines and critical fields in high 
T c superconductors, Journal of Technical Physics 40 (1999), p. 383 

[9] M. Sztyren, Bi-continuum modelling of layered structures and crystalline 
interfaces, Journal of Technical Physics 43 (2002), p. 265 



13 



